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Abstract The trace distance (TD) possesses several 
of the good properties required for a faithful distance 
measure in the quantum state space. Despite its im¬ 
portance and ubiquitous use in quantum information 
science, one of its questionable features, its possible 
non-monotonicity under taking tensor products of its 
arguments (NMuTP), has been hitherto unexplored. In 
this article we advance analytical and numerical inves¬ 
tigations of this issue considering different classes of 
states living in a discrete and finite dimensional Hilbert 
space. Our results reveal that although this property of 
TD does not shows up for pure states and for some 
particular classes of mixed states, it is present in a 
non-negligible fraction of the regarded density opera¬ 
tors. Hence, even though the percentage of quartets of 
states leading to the NMuTP drawback of TD and its 
strength decrease as the system’s dimension grows, this 
property of TD must be taken into account before using 
it as a figure of merit for distinguishing mixed quantum 
states. 

Keywords Quantum distance measures • Trace 
distance • Monotonicity under tensor products 


1 Introduction 

Quantifiers for the distance (distinguishability) between 
two density operators in the quantum state space V{1~L)- 
the space formed by positive semidefmite matrices with 
trace equal to one-are an essential and frequently used 
item in the quantum information scientist toolkit mm 

Jonas Maziero 

Departamento de Ffsica, Universidade Federal de Santa 
Maria, 97105-900, Santa Maria, RS, Brazil 
E-mail: jonas.maziero@ufsm.br 


[3j. For instance, how well a quantum system was pre¬ 
pared [5], manipulated [3] or protected [5 in an ex¬ 
periment is usually evaluated via how close (how in¬ 
distinguishable) its real state is from what one would 
ideally expect. Distance measures in T>{T~L) naturally ap¬ 
pear also in the contexts of quantum foundations |7ll8l 
0, quantum processes PHD], quantum cryptography 
mm , quantum phase transitions m, quantum speed 
limits [TMHEEHEMSD], quantum channel capacities 
[2Tj . and also in the theories of quantum entanglement 
mm, quantum discord WM®, and quantum co¬ 
herence prii2%ii!MiHn] . 

Several distance (distinguishability) measures in the 
quantum state space have been proposed in the lit¬ 
erature in the last decades. A partial list is provided 
in Ref. [31]. A few examples are the Bures’ distance, 
that is defined in terms of a similarity measure known 
as Uhlmann’s fidelity [3121155] (for a critical assessment 
regarding the use of this function in quantum infor¬ 
mation science see Ref. EH). the p —norm distance, 
with the trace distance (or 1—norm distance) and the 
Hilbert-Schmidt distance (or 2—norm distance) being 
used more frequently (see Refs. [55115611571158] and ref¬ 
erences therein), the quantum relative entropy [5511401 
mm, and the quantum Chernoff bound mm, with 
these last two distinguishability measures being defined 
operationally, respectively, in the contexts of asymmet¬ 
ric and symmetric quantum hypothesis testing. 

In this article we are interested mainly in one of 
the most popular distance measures, the trace distance, 
that is defined using the trace norm. For a Hermitian 
matrix A , the trace norm is defined and given as follows: 

\\A\\ 1 :=TrV^A = TW^='£ j \a j \, ( 1 ) 

with | u.j | being the absolute value of the real eigenvalues 
of A. We can quantify how dissimilar two density oper- 















2 


Jonas Maziero 


ators p and £ are by their trace distance (TD), which 
is defined as the trace norm of their subtraction, 

dtr{p,0 ■= Up — C||i, (2) 

and assumes values between zero and two [2]. 

This mathematical function possesses several of those 
properties required for a faithful distance (distinguisha- 
bility) measure in the quantum state space [ 2(21 : Tor 
the density operators p , £, and a, the trace distance is, 
e.g., positive semidefinite ( d tr {p , C) > 0), it is zero if and 
only if the two density operators are equal ( dt r {p , C) = 

0 O p = (), it is symmetric (d tr (p, C) = d tr {(,p)), it 
obeys the triangle inequality (dt r (p,C) < dt r {p,ct) + 
dtr (cq C)) 5 it is invariant under unitary transformations 
(dtr{p, 0 = dtr{U pU\UC,W) for UW = Id, where 
Id is the dxd identity matrix), it leads to the equal¬ 
ity dtr(p <8 a, C ® a) = dt r {p , C)i it is monotonic under 
discarding subsystems (dt r (pi, Ci) < dt r (p 12 ,( 12 ) with 
X\ = Tr 2 (xi 2 )), and it is consequently also monotonic 
under trace-preserving quantum operations (dt r {p, C) > 
dtr($(p),&(0) with 4>(x) = J2j K j xK j and = 

Id). 

Notwithstanding, it was mentioned in Ref. [44] that 
the trace distance lacks monotonicity under taken ten¬ 
sor products of its arguments. That is to say, we can 
find four density operators p , £, £, and rj such that the 
following inequalities are satisfied: 


dtr(p j C) 



(3) 


and 



,02 £02^ 


d tr (e 2 ,V® 2 )- 

(4) 


This non-monotonicity under tensor products (NMuTP) 
does not seem to be a desirable property for a distance 
measure in 'D('H). If a pair of states of a quantum sys¬ 
tem is more distinguishable than another pair of states, 
one would expect the same to hold for two identical 
and uncorrelated copies of the system prepared in those 
states. 

Two relevant questions to answer regarding this is¬ 
sue are (i) for what kind of state and (ii) how often the 
inequalities in Eqs. © and © can simultaneously hold. 
The remainder of this article will be devoted to answer 
these questions for the cases of general state vectors 
(Sec. 12.11) . for one-qubit states (Sec. 12.21) . and also for 
high-dimensional quantum systems ('Sec. 12.3|) . 

2 The non-monotonicity of trace distance 
under tensor products 

This section is dedicated to investigate such an issue 
considering some particular classes of states. Though 
we present some analytical results, much of the work 


should be numeric. We will start using general pure 
states and a two-level quantum system to address the 
question (i). In the sequence the question (ii) will be 
studied mainly with regard to its dependence with the 
system’s dimension. 

2.1 Arbitrary pure states 

Let p, £, £, and p be arbitrary state vectors on the dis¬ 
crete Hilbert space 7~L of dimension d. The trace distance 
between the pair of states x = p (£) and y = £ ( 77 ) can 
be written as [2]: 

d tr {x,y) = 2\Jl- Tr (xy). (5) 

Given that 0 < Tr (xy) < 1, dt r (p ,() > dt r (£,,p) im¬ 
plies Tr(p£) < Tr(^p). Using this inequality and the 
fact that, in the present case, 

dtr (x® 2 , y® 2 ) = 2 y/1 - [T \{xy)]\ (6) 

we see that d tr (p® 2 ,(® 2 ) > d tr ^® 2 ,v® 2 )- Thus, if all 
the states involved are pure states, the trace distance 
does not suffer from the NMuTP drawback under anal¬ 
ysis here. 

2.2 One-qubit states 
2.2.1 Collinear States 

Let us consider the special case in which the pairs of 
density operators (p, £) and (£, p) are, individually, collinear. 
That is to say, let e.g. 

p = 2 _1 (I 2 + r • cr) and ( = 2 _1 (I 2 + z • er) (7) 

with the two Bloch’s vectors being r = rn and z = ±zn, 
where h is any unit vector in R 3 and cr is the Pauli’s 
vector. One can readily show that 

dtr(p,0 = \r- F4 (8) 

For the tensor products we have 
p® 2 — C 02 = 2 -2 ((r =F z)(J -2 <8> n ■ cr + h ■ a ® I 2 ) 

+{r 2 — z 2 )n ■ cr ®fi ■ cr) (9) 

= 2" 2 ((r T z){r ± 3 + 2)P + <g> P+ 

+(r =F z ){ r ± z — 2)P_ (8 
-(r 2 - 2 2 )(P + ®P_+P_( 8 P + )), (10) 

where we used h ■ cr = P + — P_ and I 2 = P+ + P_ . It 
is straightforward applying Eq. to get 

d tr (p® 2 ,<:® 2 ) =dtr(p, C)2- x (2 + |r±*|). (11) 

We see that d tr {p® 2 ,C® 2 ) is a monotonically increas¬ 
ing function of dt r {p,Q■ Thus, for this particular set of 
states, the inequalities in Eqs. © and © cannot be 
satisfied simultaneously. 
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2.2.2 p, (, and p arbitrary states 

Even this apparently simple one-qubit case is not eas¬ 
ily tamable for analytical computations. Hence we re¬ 
course to numerical calculations via Monte Carlo (ran¬ 
dom) sampling of the quartets of states to be used. The 
computations of eigenvalues involved in this article are 
done utilizing the LAPACK subroutines (see Ref. [451 ). 
Let us start by using the Fano’s parametrization m to 
write an one-qubit density matrix x = p, £, £, 77 in the 
form: 

x = 2- 1 (h + Y^ j=1 Xj(T j ), ( 12 ) 

with x = (x-i,X2,X3), where x± = | |x ||2 sin 9 cos <f>, X2 = 
11x| 1 2 sin 9 sin </>, and x 3 = | |x| I2 cos 0. The parameters 
appearing in these equations can assume values in the 
ranges FB1I4 7] : | |x| | 2 S [0,1], 9 G [0,7r], and <j> G [0, 2w]. 
In order to obtain an uniform distribution of points 
(states) in the Bloch’s ball, each one of the quantum 
states is generated setting 

||x|I 2 = (H) 1 / 3 , 9 = arccos(—1 + 2^), <t> = 2nt 3 (13) 

with tj (j = 1, 2, 3) being a pseudorandom number with 
uniform distribution in the interval [0,1]. The Mersenne 
Twister pseudo-random number generator [35] is ap¬ 
plied to produce these numbers. By setting the Eu¬ 
clidean norm of the Bloch’s vector equal to one (zero) 
we obtain pure (maximally mixed) states. For the sake 
of illustration, the probability distribution for the val¬ 
ues of TD between pairs of randomly-generated one- 
qubit states is presented in Fig. [T] 

It is worthwhile mentioning at this point that we 
have made several tests from which we found that the 
numerical and analytical results for the TD coincide up 
to the fifteenth digit when applied to random states in 
those classes considered in Secs. 12 . 2.11 and 12.11 More 
specifically, we generated one million pairs of random 
collinear states (d = 2) and one million pairs of ran¬ 
dom pure states (see e.g. Ref. m) for each value of 
the system dimension (with d = 2, • • • , 20). The error, 
in each case, is computed by comparing the trace dis¬ 
tance obtained via diagonalization with the LAPACK 
subroutines and the value of TD obtained using its ana¬ 
lytical expression. Then the precision is established via 
the worst case error. 

We proved in the previous subsections that if the 
pairs (p,C) and (£iV) are, individually, collinear or if 
all the four states are pure, we shall have no NMuTP 
drawback of trace distance. However, as is shown in Ta¬ 
ble [I] for all the other possibilities a significant fraction 
of the one million one-qubit quartets of states randomly 
generated presented this unwanted property of TD. In 
Fig. [5] we draw an example of such a quartet of states. 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0 
Trace Distance 


Fig. 1 (color online) Probability distribution for the different 
values of trace distance (dtr(p, C) = ||r — z 1 1 2 , where £ = 
2" 1 (I 2 +z-(t)) for 10 s pairs of randomly-generated one-qubit 
states: (a) two density operators (p, £), (b) one mixed and one 
pure state (p, |£)), and (c) two state vectors (|p), |£)). The 
mean value of TD in the three cases are, respectively, 1.03, 
1.20, and 1.33. 


States generated 

Percentage 

(G) 

AG 

G max 

(p,C), (Gv) 

7.28 

0.161 

0.083 

0.475 

(p,C),(£,\v)) 

7.86 

0.186 

0.091 

0.486 

(p, IC»,(£,b?» 

3.16 

0.071 

0.038 

0.190 

(p, 0,(10, to)) 

4.06 

0.134 

0.072 

0.333 

(p,IO),(IO,to» 

3.00 

0.083 

0.044 

0.194 

(p,C),(GV2) 

8.49 

0.192 

0.098 

0.488 

(P,C), (10, Ia/2) 

20.75 

0.226 

0.096 

0.500 

(P.IOMW2) 

3.20 

0.101 

0.045 

0.248 

(p, 10), (IOJ2/2) 

7.67 

0.123 

0.037 

0.177 

(Ip), 10), (0V 2 ) 

2.63 

0.107 

0.043 

0.177 

(Ip), 10), (IOA 2 / 2 ) 

8.85 

0.169 

0.004 

0.177 


Table 1 Percentage of the 10® quartets of randomly- 
generated one-qubit states leading to the NMuTP drawback 
of trace distance. In the last three columns are presented the 
average value, standard deviation, and maximum value of the 
strength of the NMuTP drawback of TD, as defined in Eq. 
m, for each case study. 

For the sake of measuring the strength of the NMuTP 
drawback of TD , when applicable, we will define the fol¬ 
lowing quantity: 

G(p, Ci £> 9) ~ \dtr{p, C) - d tr (£,v) I (14) 

+\d tr (p® 2 X® 2 )-dtr(e\r,® 2 )\- 

The quantity G measures how far the TD is from been 
monotonic under tensor products. As G is defined only 
for those quartets of states leading to the NMuTP of 
TD, its lower bound is zero. In order to access more de¬ 
tails about the distribution of G, we shall use its mean 
value (G), standard deviation AG , and maximum value 
G max . These quantities are also shown in Table [T] A 
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Fig. 2 (color online) Example of a quartet of states for 
which the trace distance is not monotonic under taking ten¬ 
sor products of its arguments. Here r is the size of the corre¬ 
sponding Bloch’s vector, the angles are given in radians, and 
<Mp, C) = 0.80, d tr (Gv) = 0.76, d tr (p® 2 , C® 2 ) = 0.87, and 
d tr {Z® 2 ,V® 2 ) = 1-07. 



Fig. 3 (color online) Values of the strength defined in Eq. 
m for one fourth of the 72802 quartets of states leading to 
the NMuTP drawback of TD in the case study (p, £), (5, p), in 
one numerical experiment with one million random samples 
generated. The blue line is ( G) and the points in the cyan 
area are values of G in the interval [(G) — AG, (G} + AG\ (see 
Table [TJ. 


sample of the values of G for the case study (p, £), (£, rj) 
is presented in Fig. [3] 

Even with these additional informations, as can be 
seem in Table [TJ in the general case the relationship 
between the existence of the NMuTP drawback of TD 
and the classes of states involved is not an easy matter. 
For instance, starting with general states and then re¬ 
stricting one of them to be pure we pass from a percent¬ 
age of 7.28 % to 7.86 %. But then the addition of the 
same restriction for one state of the other pair reduces 
the percentage with the undesired property of TD to 


3.16 %. Several other similar nontrivial changes in the 
percentages can be identified. One striking one is that 
in the last line of the table. For four pure states there is 
no drawback, however just by putting one of the states 
in the center of the Bloch’s ball, we get a percentage 
of 8.85 %, the second higher among the classes of one- 
qubit states studied. These results stress the richness 
and complexity of the quantum state space, already for 
the composition of two of its simplest systems. Thus, 
in order to simplify the analysis, we will investigate in 
the next section the general dependence of the NMuTP 
drawback of TD with the dimension of the system. 


2.3 General one-qudit states 

In this subsection we shall study the NMuTP drawback 
of trace distance for d— level quantum systems, known 
as qudits. As there is no explicit parametrization for 
density matrices with d > 3 |Hj, we will proceed as 
follows. Let us first look for the spectral decomposition 
of a given density operator x = p, £, £, r/: 

x = H d j=i x i\ x j){ x j\- ( 15 ) 

Once the eigenvalues of x form a probability distribu¬ 
tion, i.e., 

Xj > 0 and l x j = 1 > ( 16 ) 

we can use a geometric parametrization for them m- 
Xj = sin 2 Oj^nt) cos2 0 k (17) 



3 4 5 6 7 8 9 10 11 12 13 

d 


Fig. 4 (color online) Minimum value (red bar on the left), 
mean value (green bar in the middle), and maximum value 
(blue bar on the right) of the fraction of the quartets of 
randomly-generated quantum states leading to the NMuTP 
drawback of trace distance as a function of the system’s di¬ 
mension for ten numerical experiments. We generated one 
million quartets of states in each experiment. 
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with do = 7t/ 2. The details about the numerical gener¬ 
ation of {xj} using this parametrization can be found 
in Ref. !49 . 

The basis formed by the eigenvectors of a density 
operator x , {| tcy)}, can be obtained from the compu¬ 
tational basis, {|j)}, using an unitary matrix U , i.e., 

I xj) = U\j), with j = 1, • • • , d. (18) 

There are several parametrizations for unitary matrices 
m- Here we use the Hurwitz’s parametrization with 
Euler’s angles. For details see e.g. Ref. m- 


tern’s dimension d generating one million quartets of 
states in each experiment. The mean, minimum, and 
maximum values of the percentages of the quartets of 
states leading to NMuTP of TD are shown in Fig. [4] 
In the Fig. [5] we present samples of the distribution of 
values of the drawback’s strength for some values of d. 
We can see in these figures a steady decreasing of such 
a proportion and strength as the system dimension d 
grows. 

3 Final remarks 
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Fig. 5 (color online) Samples with 5000 values of the 
strength of the NMuTP drawback of trace distance for some 
values of the system’s dimension d. The NMuTP average per¬ 
centage of the whole sample with ten sets of 10 6 quartets of 
states is shown at the side of d. The blue line indicates (G). 


d = 4; 6.51 %; <G> = 0.116 



d = 7; 5.54 %; <G> = 0.068 


d= 10; 5.15%; <G> = 0.045 





d= 13; 4.96%; <G> = 0.034 


By creating pseudo-random probability distributions 
{xj} and pseudo-random unitary matrices U, we did 
ten numerical experiments for each value of the sys- 


The trace distance has several good properties that 
rank it as one of the major distance measures between 
quantum states. Nevertheless, it also presents a poten¬ 
tial drawback, the possibility of being non-monotonic 
under taking tensor products of its arguments, that 
was shown in this article to exist for a non-negligible 
fraction of the density matrices investigated. Thus, al¬ 
though such issue seems not to be much relevant for 
high-dimensional quantum systems, it must be taken 
into account when dealing with few qubits. 

The important question that yet remains is if, in 
the cases were the NMuTP of TD is significant, it has 
some undesirable consequence for important functions 
in quantum information science. The possible implica¬ 
tions of this issue regarding, for instance, the quantifi¬ 
cation of quantum entanglement, of quantum discord, 
and of quantum coherence is an appealing topic for fur¬ 
ther researches. It would also be fruitful analyzing the 
NMuTP drawback considering other quantum distance 
measures. The obtention of a more precise operational 
and/or physical interpretation of G and its upper bound 
are also left as open problems. 
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